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Abstract. This paper revisits the theory of Y. Shikhmurzaev on forming interfaces as a continuum 
thermodynamical model for dynamic triple lines. We start with the derivation of the balances for mass, 
momentum, energy and entropy in a three-phase fluid system with full interfacial physics, including a 
brief review of the relevant transport theorems on interfaces and triple lines. Employing the entropy 
principle in the form given in m, but extended to this more general case, we arrive at the entropy 
production and perform a linear closure, except for a nonlinear closure for the sorption processes. 
Specialized to the isothermal case, we obtain a thermodynamically consistent mathematical model for 
dynamic triple lines and show that the total available energy is a strict Lyapunov function for this 
system. 


Keywords. Continuum Thermodynamics; Dynamic Contact Line; Interfacial Mass; Dynamic 
Surface Tension; Free Energy Lyapunov Functional 

2010 Mathematics Subject Classification. Primary 35Q35, 76D45; Secondary 35R35, 37N10, 
76T10, 80A99. 


1. Introduction 

The line at which three phases meet, is called a triple line; cf. Figure 1. If the phases which touch 
each other are all fluid phases, i.e. two immiscible liquids are in contact with another liquid or a gas, 
this triple line is freely deformable in space, while it is bound to move on a given surface, if one of 
the phases is a solid. In the latter case, one usually speaks about a dynamic contact line, while the 
notion of a triple line is typically used in the former setting. Both cases share many similarities and 
their modeling and analysis is closely related. In applications, wetting more often appears on a solid 
wall, i.e. the case of a contact line is more often considered. Hence, the main body of the literature 
is devoted to this case. The present paper deals with dynamics triple lines, but in such a generality, 
that analogous results are valid for the contact line situation. Nevertheless, due to the more frequent 
encounter of wetting of solid supports, the brief literature survey to follow necessarily focuses on 
contact line dynamics. 

The modeling and computation of dynamic contact lines is an active field due to the enormous 
relevance of wetting and dewetting phenomena in various technical and industrial applications; see [7], 
m and j29| for recent surveys on the field, containing also references to experimental work. Different 
modeling approaches are employed, containing in particular so-called molecular-kinetic theory (MKT; 
see, e.g., IS], m) and continuum physical theories. The latter is often subsumed under the heading 
“hydrodynamic theory” and is mostly based on sharp-interface models, while phase field models have 
also been extended to cover contact lines as in |21) . The sharp-interface hydrodynamic theory started 
essentially with the seminal paper by Huh and Scriven [20] in which the fundamental problem of the 
inconsistency between a moving contact line and a no-slip condition at the fixed wall has been analyzed 
and shown to lead to a non-integrable stress singularity; cf. also [16] and, for a more rigorous mathe¬ 
matical treatment, [25] and [30]. Consequently, subsequent models always rely on some “relaxation” 


Date: April 21, 2015. 


1 



2 


D. BOTHE AND J. PRUSS 



Figure 1. Configuration of the phases and interfaces at the contact line. 

at the contact line, and the most common way to remove the stress singularity, as already proposed 
in [20j . is to introduce Navier-type slip close to the contact line. Besides this complication, the main 
extension of the standard two-phase Navier-Stokes system consists of the prescription of the dynamic 
contact angle Od, i.e. the angle which is formed between the fluid interface and the solid support, as 
a function of the contact line speed. At this point it is to be noted that the contact angle changes 
its value under dynamic conditions, while the equilibrium contact angle 9^ is usually assumed to be 
governed by Young’s law, i.e. 

cos{9e) = (1) 

in case of a liquid wetting a solid surrounded by a gas phase, where the superscripts stand for gas (g), 
liquid (1) and solid (s) and a denotes the interfacial tension of the respective interface. Based on the 
classical experimental studies in m, the general form of the relation between and the contact line 
speed is given by the heuristic relation 

9d = /Hoff(Ca + /Hoff(^e))’ (2) 

where Ca denotes the Capillary number given as Ca = rjU/a^^ with rj the dynamic liquid viscosity and 
U the contact line speed. Several concrete correlations have been established for different materials 
and certain wetting scenarios like “Tanner’s law” [3T] or the correlation of Jiang et al. [22]. Theoretical 
investigations using the hydrodynamic theory identified three length-scales near the contact line: an 
inner region in which the fluid interface is essentially planar and touches the solid support at the 
equilibrium angle; a mesoscopic region in which a significant bending of the interface can occur; an 
outer (macroscopic) region in which the contact angle attains a different value, the so-called apparent 
contact angle. The hydrodynamic theory provides relations for the dependence of the contact angle on 
the distance from the contact line especially in the mesoscopic region; see m, m, m Knowledge 
of this dependence is very useful for numerical purpose, both as a subgrid-scale model to reduce the 
necessary resolution at the contact-line and in order to neutralize the inherent mesh dependence of 
numerical solutions due to the typical under-resolution of the smallest length scales in the contact line 
region; for the latter, see |T| and m- 

While the hydrodynamic model can describe many wetting processes at least qualitatively, in par¬ 
ticular concerning the observed dynamical shapes of attached droplets moving on a wall, say, it does 
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not capture the full physics of a dynamic contact line. One important deviation is the internal flow 
field in the wetting liquid close to an advancing contact line, which experimentally is known to be a 
rolling motion (uni, m), but is a sliding motion in numerical simulations using the above model. 
Moreover, there is experimental evidence that the relation between the dynamic contact angle and 
the contact line speed is more complicated and not of such a simple local nature; cf. [8]. For further 
discrepancies between experimental observations and the hydrodynamic model see m- A very inter¬ 
esting approach to overcome these short-comings has been introduced by Y. Shikhmurzaev in |26) : see 
also m- The approach there also employs continuum physics, but accounts for the aspect of interface 
formation and disappearance at the contact line. A crucial point for the model development then is 
to include enough interfacial thermodynamics to allow for a non-constant interfacial tension, governed 
by a surface equations of state on all involved interfaces. For this purpose, the mass contained in the 
interfacial layer has to be balanced separately, since it encounters different forces compared to within 
the bulk phases and it is this mass density which determines the surface pressure, i.e. the surface 
tension. In the considered sharp interface/sharp contact line model, the interfacial mass is lumped 
into an area-specific mass density and the model is extended to cover the evolution of this interfacial 
mass density by appropriate balance equations on the moving surfaces. This model has proven a great 
potential to explain several physical phenomena like wetting, coalescence, cusp formation and the 
break-up of liquid threads; cf. [27], m and the introduction in |6]. The interface formation model of 
Shikhmurzaev has been based on the continuum thermodynamics of fluid interfaces developed in |5], 
|3] , but with the sensible aim to formulate the most simple model which is able to describe the wetting 
process with dynamic contact angle and rolling motion close to the contact line with the material 
properties modeled via bulk and surface free energies but without a heuristic relation between contact 
angle and contact line speed. Several years after the fundamental paper |26j appeared, Billingham 
in |B] pointed out that one further condition at the contact line has to be added, and he employed a 
condition provided by Bedeaux in |1] which relates the rate of mass transfer from one surface into the 
other to the difference of the surface chemical potentials. We will come back to this point in the final 
remarks at the end of this paper. 

In the brief survey above, the topic called “contact angle hysteresis”, referring to the appearance 
of a full interval of possible contact angles in the static case which spans the range from the angles 
observed for (infinitely slowly) advancing and receding contact lines, has not been touched. This 
phenomenon seems to be similar to dry friction between solids and, in fact, the notion of contact line 
friction is also present in the literature on the molecular kinetic theory of contact lines. For this topic, 
we refer to [TO], [29| and the references given there. 


2. Integral Balances 

We consider a region G C filled with three bulk phases {k = 1,2,3), separated by interfaces 

S^(t) {k = 1,2,3) which meet at a common triple line C. As an example, imagine a liquid phase 12^ 
in the form of a water droplet sitting on another liquid, say oil, which forms bulk phase and being 
surrounded by phase composed of air. Then, for instance, denotes the oil-water interface, Y? 
the interface between the air and the oil and the air-water interface. The deformable and free 
bounding curve at which all three interfaces meet is the so-called triple line C. As a related but 
somewhat different case, consider again a liquid phase in the form of a droplet, but now sitting on 
a solid support, which forms bulk phase and being surrounded by gas phase Vfi. Then two out of 
the three interfaces are fixed and the triple line is the set of all points where the gas-liquid interface 
meets the solid support. In this case, one usually calls C the contact line which now has reduced 
degrees of freedom due to the solid support. We focus on three-phase fluid systems with a common 
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triple line and assume that the interfaces meet at angles different from 0 and vr; we shall refer to this 
as the non-degenerate case. 

We start with the integral balance of a generic extensive quantity which is present in the bulk phases 
with specific density cj), on the interfaces with specific density and on the triple line with specific 
density (jf. Hence pcf), p^4>^ and p^cjf, respectively, are the volume-, area- and line-specific densities, 
where p, and p^ are the mass densities. If a specific bulk phase or interface is considered, we write 
Pk(t>k or Pk(t>k^ respectively, for the respective density. With this notation, the generic integral balance 
for a fixed control volume V G G reads as 


A 

dt 



p(j)dx -|- 



do -|- 



/ {P<pv + j) -ndo- / {p^<t> 

IdUv JdT,v 


+ j^) ■ Ndl - 


+ j^) ■ vdP 


' dCy 


(3) 


+ [ fdx+ [ f^do+ [ fdl. 

J'J ^Y JCy 

In Q we let dx, do and dl denote the volume, area and line measure, respectively. Moreover, dP 
denotes the point (i.e., counting) measure. 

Here as well as throughout the paper we use the following condensed notation. First, 

3 3 

Hy := IJ Sy := [j with := n V, := n V, 

k=l k=l 


which are all time-dependent sets. We assume a single triple line, hence Cy := C n V. The detailed 
version of Q then reads as 

3 ^ 3 

Pk(Pkdo+ I p^4>^dl 
ICv 


A 

dt 


■ 3 „ ^ r 

V / Pk(t>kdx + V / 

^^ /ofe ^ ^ lyk 


3 » r r 

V / {Pk(t>kVk+jk)-n’^do-y2 {Pk(l^kVk+jk)-d^^dl- {p^(lfv^+f)-vdP 

tZi da^v kP^ da^’^ Jacv 


k=ida^^y 


On On n 

tV / / fkdo+ fdl. 

* ^ I C)k * ^ /v'fe ]p 

k=\'^^v 

For better readability, we use the condensed notation whenever this is reasonable. 

We apply this balancing to the extensive quantities mass, momentum, energy and entropy. The 
corresponding integral balances read as 

mass balance. 


dt 



pdx + 




/ pv ■ ndo 
ddflv 



p^v^ 


Ndl 



• vdP. 


(4) 


momentum balance. 
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f pv{v-n)do— f p^v^{v^ ■ N) dl — f p^v^{v^-v)dP 
Jonv JdT.v JdCv 


(5) 


+ [ S-ndo+ [ S^-Ndl+ [ S^-udP+ [ pbdx+ [ pH^do+ [ p^b^^ 
J d^\r J dYj\A J dCv J ^\r J C\r 


(fb^ dl. 

^ dCy JQy '^'^y JCy 

Here S, and denote the stress tensor in the bulk phases, on the interfaces and on the triple 
line, respectively, and b, b^ and b^ are the specific body forces. 


dt 


energy balance. 

[ Pie + Y)dx+ f p^{e^ + ^^'^)do+ j p^(e^ + -^^)c 

_Jq,y ^ Jy^y ^ J Cy ^ 

- [ p{e + Y'jy . n da - f p^{e^ + ^ )v^ ■ N dl - [ p^{e^ + ^ dP (6) 

Janv ^ JdT,v ^ JdCv ^ 

+ f {v-S — q)-ndo+ f {v^ ■ — q^) ■ N dl + f {v^ ■ — q^) ■ v dP 

Janv JaT.v Jacv 


7 S 


+ / pv-bdx+ I p^v^-b^do+ I p^v^-b^dl. 
J Q,y 'J Yy Cy 


Here q, q^ and q^ denote the heat flux in the bulk phases, on the interfaces and on the triple line, 
respectively. Note that energy sources due to radiation have been omitted in Q. 

entropy balance. 


A. 

dt 


I psdx+ I p^s^doP I p'"s'" dl 

' ^Y 'J Yy Cy 


^cx, 


f {psv + ^) ■ ndo — [ {p^s^v^ + d>^) ■ N dl — f {p'"s'"v'" + <1>^) • i' dP 

Janv JdT,v JdCv 


(7) 


+ / Cdx+ ydoP Y dl. 

J ^Y J Yy J Cy 

Here <h, and denote the entropy flux in the bulk phases, on the interfaces and on the triple 
line, respectively, while C, and C!" are the corresponding entropy productions. 

Remark. Note that the internal energy density as well as the entropy density can be positive even 
if the area- or line-specific mass densities are considered to be zero. In other words, in the limit as 
p^ -a- 0-|- or -a- 0-I-, products such as or p^s^ may converge to strictly positive limit densities, 

i.e. 


s s 
p e 






c c c c 
vr , p^ d" ^ 


with non-vanishing densities has to be allowed for. Otherwise, for instance, the surface 

tension for a fluid interface with zero surface mass density would automatically vanish. 

3. Transport Theorems 

The derivation of local versions of the balance equations follow by application of appropriate trans¬ 
port theorems and subsequent localization. The following transport theorems will be employed. 
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Volume transport. In the general setting described above, let V C be a fixed control volume in 
G, let S be short for with the time-dependent interfaces and ns = nj^k the unit normal 

field on with an arbitrary fixed orientation. Let denote the speed of normal displacement 

of The latter is a purely kinematic quantity, but it is related to the barycentric velocity of the 

interfacial mass via Vs = • ns- Moreover, given any bulk field (f), the jump of at S is defined by 

the jump bracket [-J according to 


10] (t, x) := (0(L X + hn^) - x - hns)) ■ 


( 8 ) 


With these notations and for the specific control volumes mentioned as well as for sufficiently smooth 
fields, it holds that 

^ [ 4>dx= f dt(j)dx - [ 101 Vs do, (9) 

dt Jv Jv\j: Jt.v 

where T^vit) := S(t) n V. 

Surface transport. In the general setting described above, let V C be a fixed control volume in 
G. Then, for sufficiently smooth fields, it holds that 

d 


dt 


(jp do= [ - 0^nsVs) do+ [ cj)^ Vgsv 

JdT.v 


dl. 


( 10 ) 


Here denotes the time derivative along a path that follows the normal motion of S(-) and ns := 
divs(—ns) is twice the mean curvature. Furthermore, VdT,v is the normal (relative to the boundary 
of Ey) speed of displacement of (9Ey(-) (in the plane tangential to S). 

Let us note in passing that the derivation of the local balance equations can be done with special 
control volumes such that the outer normal ny satisfies ny T ns on dV n E. For such control volumes 


the boundary contribution, i.e. the final term in (10), vanishes. 


Line transport. For sufficiently smooth fields, it holds that 

(lfdl=f + (l)‘^dwcv‘^) dl + [ cf {VdCy - v‘^ ■ u) dP. (11) 

dt JCv JCv ^t JqCv 

Here ^ denotes the Lagrangian derivative, following the triple line along a path with velocity and 
Vqcv is the normal (relative to the end points of Cy) speed of displacement of dCy- Recall that u is 
the outer normal to the curve Cy in its end points (cf. Figure 1) and that dP denotes the point (i.e., 
counting) measure. 

Remarks. 1. The transport theorems above appear rather different. Actually, they can all be brought 
into the same form as the line transport theorem. In case of surface transport, this follows directly 
from the relation 


Dt 


= 9* 0 V ■ Vs0 


for the surface Lagrangian derivative. Note that div^u^ = divsu0 — ksVs, hence 


dfctP - 






u0) -b 0^div 


and then, by the surface divergence theorem, equation (|10|) implies 

(p^ ao = I ( 

IJ^v 


dt 


(fP do= [ jf. + 0^div do+ [ cjP (VaSv - ■ N) dl. 


( 12 ) 
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To bring the volume transport formula ([^ into the same form, one first observes that ([^ combines the 
transport formulas for both bulk phases which meet at the considered interface. If two bulk phases 
are separated by an interface S(t), then a simple variant of the Reynolds transport theorem 

yields 

d 




D(j) f 

( 7 ^ + (pdwv)dx+ 6 {Vq^± - vn) do, 
Jdnt ^ 


(13) 


where n is the outer unit normal to fly; note that the latter coincides with ±ns on Sy. 


2. An equivalent form of (12) reads as 




/ 


Dt 


+ (fPdiv Y.v^) do — / 


■ nv 


/asv 


a/ 1 - (ns • nvY 


dl. 


(14) 


where ny is the outer unit normal to V. For this purpose, one first uses elementary geometry to 


compute VdEy = -Vs 




/l-(ns-ny)2 


Since -/l — (ns ■ ny)^ = N ■ ny and {N,ns,T} with r a unit 


vector tangential to clSy (hence also to dV) is a local orthonormal basis, the equation (14) follows 
from 


JV - VffSy = 


1 


-/l - (ns ■ nyy 


■ N){N ■ ny) + {v^ ■ ns) (ns • ny)^. 


The relation from (12) has been given in |12) . while the variant (14) can be found in Chapter 3 in 
see also the appendix in [ 2 ]. 

3. Below we will also use variants of the above transport theorems with built-in mass balance. 
These read as 

d f f ^ 

— / p<pdx= / p-—dx + 


dt Jy JV\T, 

with := p^{v^ — v^) ■ ns on S, and 


Dt 


do — p(j)v ■ n do 
Sv JdV 


(15) 


A 

dt 


p^(j)^ do = 


iv 






- do+ p^4P (Fasy -v^ ■ N) dl. (16) 

JOHy 


Dt 

4. Local Balances 


Application of the transport theorems and localization yields the following local balance equations 


bulk phase. 

dtp + div (pv) = 0, (17) 

dt{pv) + div {pv iSi V — S) = pb, (18) 

dt{pe) + div {pev + q) = S : Vv, (19) 

9/ps)div (psu-h 4)) = C- ( 20 ) 

These are the well-known balance equations in a bulk phase, 
interface. 

dfp^ + divsip^v^) -h Ip{v - v^) ■ ns] = 0, (21) 

df{p^v^) + divs(/0^f^ - S^) -h l[pv ® {v - v^) - S) ■ ns| = P^h^, (22) 


{ 5 ~! \ 9 ! 

df{p^e^) + diYs{p^e^v^ + q^) -h l(^p{e + ^ ^ ){v - v^) - (v - v^) ■ S + q^ ■ nsj 
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= (23) 

df{p^s^) + divs(/9^s^z;^ + + |(ps(t> - v^) + d>) • nsl = C^- (24) 

Observe that the jump terms always appear with as a factor. Therefore, these terms are invariant 
under re-orientation of the interfaces. Actually, the notion of a “jump condition” for these terms can 
be misleading. Note that |/ • — f~ ■ n~ if the interface separates two bulk phases 

with outer unit normals n^. Hence, if / denotes a bulk flux, the term — |/ • ns] describes the total 
rate of transfer from the bulk phases to the interfaces due to these fluxes. For the derivation of closure 
rates below, the explicit form of this term is to be used since two binary products are involved. 

triple line. 

/ + divc(//) + |p^(u^ - /) • iVl = 0, (25) 

af (/u^) + divc(// ® - S^) + |(/9^u^ ® {v^ - v^) - S^) ■ = p^b^, (26) 

^(/e^) + divc(/e^u^ + /) -b l(/o^(e^ + - )(u^ - v^) - {v^ - v^) ■ • iV| 

= : Vc/, (27) 

(/s^) + divc(//u^ + ^^) + |(p^s^(^^^ - v^) + $^) • iV| = C^. (28) 

Here the triple bracket l-J is dehned exclusively for quantities of the form f'^-Nhy means of 

3 

lf^-Nl = -Y,fk-N^ onC, (29) 

k=l 

where the sum runs over all interfaces which meet at the triple line and := /|sfe- Let us briefly 
explain the appearance of such terms, e.g., for the mass balance (25). The transport relation (16) for 
(j)^ = I yields the boundary contribution of the interfacial mass balance as 

3 


[ {Vaj:, -v^-N)dl = ^ [ pI - vf • AT^) dl. 

JdT.v ^ 

The boundary of TJy is (S*^ n dV) U Cy, hence 

[ p^{V9^,-v^-N)dl = Y, [ pUV9j:^^-Vk-N^)dl+ [ ^ (/- • iV^ d/. 

rr J&nav ^ Jcv rr 


I9T.V 

Employing the condensed notation, this becomes 


/asv 


{Va^v -v^-N)dl= [ P^ (Hasv - ■ N) dl - [ |p^ (/ - v^) ■ iV] dl. 

dsnSP JCv 


5. Entropy Production and Closure Relations 


The entropy principle states that every admissible closure for the entropy flux is such that the re¬ 
maining entropy production is a sum, running over all dissipative mechanisms, of binary products. The 
entropy production is non-negative for any thermodynamic process, i.e. the entropy inequality holds. 
The system is in equilibrium, if and only if the entropy production vanishes. For more information 
about the employed entropy principle see m- We are going to apply this for bulk, interface and triple 
line in a fully analogous manner; the details will only be explained for the bulk case. We consider the 
simplest class of bulk, interface and contact line materials for which the entropy density is assumed to 
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be a concave function of temperature and mass density, only. We hence employ constitutive relations 
of the form 

ps = h{pe,p), p^s^ = h^{p^e^,p^), p^= hP{p^, p^) (30) 

with concave functions h, hP and hP. We furthermore define the (absolute) temperature and the 
chemical potential in the respective phase as 


and 


1 dh 

1 dh^ 

1 

dh^ 

(31) 

T d{pe) ’ 

^(^SgS)’ 

pc 

d{p^e^) 

p dh 

p^ dh^ 

~ dp^ ’ 

/ 

pc 

dhp 
dp^ ' 

(32) 


We insert the constitutive relation (poj) for the entropy density into the respective entropy balance, 


use the chain rule employing the definitions (31) and (32) and eliminate all partial time derivatives 
by means of the other balance equations. The resulting terms are grouped in such a way that only 
a single full divergence appears, which contains in particular the entropy flux, all terms with the 
velocity divergence as a factor are collected and all remaining terms are grouped to form a sum of 
binary products. 

bulk phase. The procedure above yields 


( = div ^{pe + P - psT - pp) dSNv + q - vl- + : Vn, 


(33) 


where P := — ^trS* is the mechanical pressure and S° := S + PI, with I denoting the identity tensor, 
is the traceless part of S. We will assume throughout this paper that the material in all phases does 
not support local densities for angular momentum (so-called couples). Hence the balance for angular 
momentum implies that all stress tensors which appear are symmetric; note that all stress tensors are 
formulated in the embedding three-dimensional Euclidean space, i.e. are symmetric 3 x 3-tensors. 

Evidently, the simplest closure for the entropy flux in order to fulfill the entropy principle is := ^, 
which is the standard choice for single component materials. This leads to the reduced entropy 
production, being the desired sum of binary products. Exploiting the symmetry of S, we obtain 


C = -^(pe + P - psT - pp) dWv + q ■ V^ + ^5° : D°, 


(34) 


where D := |(Vn -|- (Vn)''') is the symmetric part of the velocity gradient and D° its traceless part. 
The dissipative mechanisms associated with these binary products are “volume variations”, “heat 


conduction” and “viscous shear”, in the order of their appearance in (34). The simplest linear (in the 
co-factor) closure without cross-effects leads to the relations 


pe + P — psT — pp = —A div v 

with A > 0, 

(35) 

^ 1 
q = aV- 

with a > 0, 

(36) 

S° = 2r]D° 

with 1 ] >0. 

(37) 


Note that the closure parameters X,a,r] are allowed to depend on the basic variables, say {p,T). 
Hence, in particular, the heat flux closure is equivalent to Fourier’s law. For consistency with standard 


notation, we use 2ri instead of rj, above. At this point, an explanation concerning (35) is at order: The 
only quantity which requires a closure is P = — ^tr5. For a stagnant fluid, equation (35) reduces to 
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pe+p—psT—pjjL = 0, where p denotes the pressure at equilibrium. We therefore let the thermodynamic 
pressure p be defined by the Gibbs-Duhem relation, i.e. by 


pe + p — psT = pp. 


Then P = p + tt with the non-equilibrium pressure contribution vr and (35) becomes 

TT = —A div V. 


(38) 


(39) 


The irreversible pressure contribution vr is due to volume variations and the linear closure to model it 
reads as tt = —A divu. Let us note that the thermodynamic pressure p from (38) satisfies the Maxwell 
relation P = P^^ with the free energy if} = 'ijj{T,p) := e — sT. Alternatively, one can define p by 
the latter relation and obtain the Gibbs-Duhem relation (38) as a consequence. Note also that the 
entropy production (34) can now be written more concisely as 

: D, (40) 


C = g-vi + i5“ 


where the irreversible stress part is defined as = —ttI + S°, but it is important to notice that the 
last term represents two independent binary products. 

interface. The same line of arguments leads to 

..s 

J S S I S S S S /' /11A 

and p e +p —psT =p p (41) 






as well as 




^ 1 1 s V s 

Vs^-^vr divsu 


+ : D 


S,o 


+ 


- u^)|| • (5 • nE)||] + [(^ - (m(e + -) + q-n^\ 


(42) 


l \ rp rpY, ' rpY, V 


T 




1 Av — v^) 


E\2 






riY. 


ns) m 


recall that m = p{v — v^) ■ n-£. Here vr^ is the irreversible part of the interface pressure defined 
via + p^ = — ^tr5^ with the thermodynamic interface pressure p^ from (41)2- Moreover, = 
^Is(Vsu^ + (Vsu^)''')/s is the symmetric interface velocity gradient, D^’° its traceless part and 
Iy, = I — riY, ® riY, denotes the surface projector, also called surface identity. 

The dissipative processes associated with the binary products in (42) are, in the order of their 
appearance, interfacial heat conduction, area variation, interfacial shear, one-sided slip between the 
interface and a bulk phase, heat transfer to and from the interface and, finally, mass transfer to and 
from the interface. The following closure relations result by assuming linear relations between the 
corresponding co-factors with one exception: the mass transfer to or from the interface, i.e. the ad- 
and desorption processes m = will be modeled using a non-linear relationship in analogy 

to the modeling of chemical reactions; cf. m- 




s \S 1 - s 

TT = —A dlV yV 


with > 0, 


with A^ > 0, 


(43) 

(44) 


= 2r/^D^’° 

with > 0, 

(45) 

lA^{v - u^)|| (571^)11 = 0 

with 0^ > 0, 

(46) 
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^ - 7 ^+^'^{p{e+-^){v-v^) + q) - ns = 0 


or In 


m 


fi 1 Av-v^)^ , 




with 5^ > 0, 


with > 0. 


(47) 


(48) 


The closure relation (48) employs the decomposition rh = Note that (48) only fixes the ratio 

of ad- and desorption, while one of the rates needs to be modeled based on experimental knowledge or a 
micro-theory. The simplest choice is to assume a desorption rate according to with > 

0. Observe also that (48) is an implicit equation regarding m, since = (u^ — u^)|| -|- mA/. 

At this point it should be noted that the closure relations above are given in a condensed notation: 
relations (43), (44) and (45) are employed for every interface TA [k = 1, 2,3) separately with respective 
transport coefficients, while the transmission relations (46), (47) and (|48|) apply to each interface in 


combination with any of the two adjacent bulk phases. In total, the closure relations hence yield nine 
conditions at each of the three interfaces. 

triple line. Since the triple line is one-dimensional, the contact line stress tensor satisfies = —P^Ic 
with the mechanical line pressure := — trS'*' and the line projector defined by Ipw = {w,t)t with 
r a unit tangent field on C. By the same procedure as above, we obtain the following identities, where 
P^ = + vr*', and fnp = p^{v^ — v^) ■ N, i.e. rri^ = p^{vf — v^) ■ for A: = 1, 2, 3. 


as well as 


^ y and 

T'~- 


(f = q^ ■ - 7^71^ dive/ + ^|(/ - /)||| • {S^ ■ iV)|||l 


(49) 




1 


j'C j'C 
1 


s 


TC 

,c 






■m 


(50) 


/I _ /I I J_ 

j'S J'C ' jC 




m 


2 p- 

Above, the notation (•)||| denotes the component tangential to the triple line and := —vr^/s -|- 

S^’°. In analogy to the interface we obtain the following closure relations for the dissipative processes 
on the triple line. 

/ = /Vc^ 


TT^ = -X^ divcm 


/(u^ — / 
1 


1 

' 

/ In 


+ (5^iV)|||=0 




m 


jC 

E.ad 




C\2 


m 


S,de 


pP 1 Av^ — V^) 
~ 2 


-N ■ 



with / > 0, 

(51) 


with > 0, 

(52) 


with / > 0, 

(53) 

r = 0 

with / > 0, 

(54) 

cS.irr 

/ .N) 

P 

with / > 0. 

(55) 

= p^{v^ — /) 

• N as fnp = fnP'^'^ — 

• y Hfs . 

m ’ IS 


employed. The relation (55) governs the ratio of ad- and desorption at the triple line, while one of the 
rates needs to be modeled based on experimental knowledge or a micro-theory. Below we assume the 
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desorption rate to be given by with > 0. The transfer relations (53), (54) and 


(55) exist for every combination of the triple line with one of the interfaces, of course with individual 


transfer coefficients. 

To complete the model it remains to fix free energy functions for the bulk phases, the interfaces 
and the triple line. This will only be done for a reduced model below. 

6. Isothermal Case with Vanishing Triple Line Mass 

We consider the limiting case of isothermal conditions, i.e. the internal energy balances are replaced 
by a known constant temperature field; in particular, we have T]s = and = T^. We also reduce 
the model complexity by neglecting the mass and inertia on the triple line. Moreover, we neglect 
any irreversible stress contributions both on the interfaces and on the triple line. For consistency 
with the notation in interfacial science, we do not employ the surface and line pressure, but rather 
let = 7 ^/s and = 7 ^ic with the interface tensions 7 ^ and the line tension 7 ^. Observe that 
this means = 0. Because of zero triple line mass and isothermal conditions, we assume the line 

tension to be constant. 


bulk phase. 

dtp + div {pv) = 0 , 

dt{pv) + div {pv ®v) = div S + pb, 


(56) 

(57) 


where the stress is given by S' = {—p + Xdivv)I + 2r]D° according to (37) and (39). In the compressible 
case, an equation of state in the form p = p{p) (with a strictly increasing function p(-)) is to be added 
according to the specific fluid under consideration. 

interface. We again use the abbreviation m = p(v — v^) ■ njj. Then 

dpp^ + divs(/9^u^) + [m] = 0, (58) 


dtiP^v'^) + divs(p^u^ (g) v^) + |um] = [S' • nsl + div^S^ + p^h^. 






Note that 


div + Ve 7 ^ 


in the considered case without surface viscosities. 

This is complemented by the constitutive transmission conditions 


/3 (v-v )|| + (S'nE)|| =0, 

ad 


El rn 
a In 




= p-p + 


gzrr 

ns-ns, 


(59) 

(60) 

(61) 

(62) 


2 p 

where > 0. In addition, the material dependent interface free energy function is required. The 

latter determines especially the interfacial equation of state 7 ^ = 7^(/9^) and we assume that 7 ^ is 
a strictly decreasing function (i.e., the interface pressure depends strictly increasing on the interface 
mass density). 

triple line. In analogy with the interface-related notation, we use as before the abbreviation 

:= p^{v^ — u^) • N, 

i.e. = pf{v^ — v^) ■ for A: = 1, 2, 3. Due to p^ = 0, the triple line mass and momentum balances 
become 


|rh^| = 0, 


(63) 
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(64) 


= +7‘^divc/c- 

This is complemented by the constitutive transmission conditions 

^Riii = ^'7111 = ^3;iii =^4’ 

- / + {{vf - /) • N^f/2 = 0 {k = 1,2,3). (66) 

A few comments are at order: for simplicity, we consider the no-slip condition (u^ —u^)||| = 0, but note 
that the barycentric triple line velocity is undefined for a triple line with zero mass. We consider 
as the kinematic velocity of the contact line C. The chemical potential is determined by one of 


the equations in (66), so actually only two equations remain there. Also, observe that —v^ T r, ns, 
hence 

— v'"\N) N, (67) 

as C{t) C Sfc(t) for all times. 


7. Thermodynamical Consistency and Equilibria 

For this reduced isothermal model we show that the total available energy, i.e. the sum of the total 
kinetic energy and the total free energy is a strict Lyapunov function in case of vanishing body forces. 
We hence let 


+ Ip ) do+ dl, 


( 68 ) 


Ea(t) = f p{Y + 'iP)dx+ j p^{- 
where G is the total domain. We are going to show 

Theorem 1. Let {p,v, p^,v'^,T,,C) he a classical solution of the model from i.e. a classical 

solution to (56), (57) with S = {—p + Xdivv)I + 2pD°, where p{p) is strictly increasing in p, X,t] > 0 
and 6 = 0, (58), (59) with = 7^7s, where 7^(p^) > 0 is strictly decreasing in p^, and 6^ = 0, ([6l|) 


with > 0, (62) with > 0, (63), (64) with 7*' a positive constant, (65) and (66). We also assume 


that this solution is non-degenerate at the contact line, i.e. the interfaces meet at angles different from 
0 or TT. At the outer boundary, we assume v ■ n = D, v ■ Sn = 0 on dG and ■ N = Q on S n dG. 


Then the total available energy Eq from (68) is a strict Lyapunov function. 


Proof, (i) Let {p,v, p^ ,v^ ,T,,C) be a classical solution of the model from Section]^ For the 
bulk contribution, we first apply the transport relation (15) and use the momentum balance ([T^ 
to eliminate p%i. We then exploit ip = ip{p) with ip'{p) = pj and use the mass balance ( flT] ) to 
eliminate Application of the two-phase divergence theorem for partial integration in the form 


V ■ divS dx = 


IG 


f v-Sndo— f SiVvdx— f |n • S'ns] do 
JdG Jg Jt. 


yields 


4 

dt 


'G 


f 

p{gp + i’) dx = - I p{—-\-ip)v ■ ndo-\- 
^ JdG ^ 


V ■ Sndo 


’dG 


- / 

Jg 


:’Vvdx— / |n • 5ns] do-|- 

JS 


|m(V’-E —)] do. 


(69) 


For the interface contribution, we first apply the transport relation (16) and use the momentum 
balance in the non-conservative form 

^ -E |(u - u^)rhl = [5 • ns| + divs5^, 


Dt 
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which follows from (22) and (21), to eliminate . Next, we apply the surface divergence theorem 

for partial integration of •divs5'^, employ with {p^) = /{p^Y = — 

(21) and the interface Gibbs-Duhem relation (|41|)2 to obtain 

4 [ do = [ +tp^){Va^-v^ ■ N)dl+ [ v^S^-Ndl 

Ot Jy. ^ JdT, ^ JdT. 


- / 5^’” : 


do + f i;^-|5nsldo+ f - ■ v — p^)rhldo. 

JE is 2 


(70) 


For the triple line contribution, we apply the transport relation gi, For constant 7 ^ this yields 

^ y 7*^ d/ = J 7‘'divci;‘' dl = j 7^/c : Vcu*' dl = —7*^ J ■ div cic dl. 

Employing the identities ( [ 6 ^ , (70) and ( [tT] ), we obtain 

Ea=/ + ■N)dl+ [ v^S^-Ndl 

JdT. 2 ias 

- [ s^'^^-.Ddx-J do - J |(i;-u^)|| • ( 5 ‘’'''ns)||ldo 


(71) 


IG 


r ( 1 ) — 7;^')2 cirr p 

+ / lip- + - -z- riY ■ -• ns) rh] do - 7 ^ / • divcic dl. 

JT. P Jc 


(72) 


Inserting the constitutive relation S'^ = 7 ^/s and exploiting the assumptions u • n = 0, u • Sn = 0 
dG and = 0 on S n dG, we get 

Ea = - f S^'^^-.Ddx- f |(u - u^)|| • (5‘’''’nE)||] do 


on 


>G 


+ / lip — P^ + 


iv - S'*"" 

---ns •-• ns) nij do 

z p 


,iv^f 


+ 


+ V’^)ni^l dl — f | 7 ^u^ • A^l dl — ■ div^/c dl. 


(73) 


Expanding | 7 ^u^ • A^| as | 7 ^(u^ — + v^) ■ A^|, exploitation of (64) allows to rewrite the triple line 

contribution as 

,,Sl2 


y dl- 


Using (u‘')^|rh^| = 0 due to (63), where i;|j| is given as the well-defined tangential part by (65), 
we see that ( [TS] ) implies 

Ea = - / S*"" :Ddx- [liv- u^)|| • (S‘""ns)||l do 
Jg ds 

" s , iv-v^)^ S*"" 


+ J Up — P^ d- 


riY ■ -• ns) ni] do 

P 


E , iv^-V^f 


l(/^^ + 




2 


)rh^| dl. 


(74) 
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To come to the final representation of E^, we have to write out the jump brackets !•] and 


We 


start with the triple line contribution and have, by (63), 

3 






= - 


- / + 




k=l 


Since — — {{v^ -v^)- by (67), the triple line contribution vanishes due to the constitutive 

assumption ( 66 ). Insertion of the other constitutive relations, i.e. (37), (39), (61) and (62), finally 
leads to 


Ea = — / X{divv)'^ dx — / 2riD° : D° dx 

Jg Jg 

- j -v^)‘^^do- j do 


- / a 


' (log — log do 


- / a 


s,- 


(logm-’'^^ - logm-’^'^) do. 


(75) 


Notice that, according to our condensed notation, the integrals over S are to be taken over the three 
interfaces and the notation (•)^ then denotes the respective one-sided bulk limits. Evidently, (75) 
shows that E^ is decreasing along classical solutions, i.e. Ea is a Lyapunov function. 

Next we want to characterize the equilibria of the problem, proving at the same time that the total 
available energy Ea is a strict Lyapunov functional for the system. To this end assume that we have 
a solution where Ea is not strictly decreasing at all times. Then there is an interval J = (ti,t 2 ) where 
Ea is constant, hence dEa/dt = 0 in J. This implies, by (75), 


divu = 0, D° = 0, 


V\\ = '^11 = 'I’ll 


771 + = m = 0 , 


as A, r?,> 0 by assumption. This yields D = 0, as well as |u] = 0 on S, which by Lemma 
1.2.1 of the monograph |24) implies u = = 0. Next, investigating the equations for the bulk, we 

see that dtp = 0 and Vp = 0 , which implies that p is constant in the phases, as pk is by assumption a 
strictly increasing function of pk- 

In the next step, we look at the equations on the interfaces. By the definition of we obtain 


0 = 


(u* — v^) ■ 


s s 
• n , 


hence ■ vP = 0 which yields = 0. Then the mass balance on S implies cl|"ps = 0 on J. 


Furthermore, v = 0 and p constant yield constant, hence = pP is constant by (62). This shows 
that p^ is constant, as pP is strictly increasing with p^. To see the latter, recall that p^ = +p^/p^ 

and, hence, (p^)^(p^) = {p^)'{p^)/p^ > 0. This shows further that 7 ^ is constant. Looking at the 
stress transmission condition this further yields kP constant on each of the surfaces more precisely 
we obtain = |p], i.e. the Young-Laplace law holds on each of the surfaces 

In the final step, we consider the equations on the contact line. Here we have U|j| = 0 by (65), as 
well as 


= {v^\N’^)N^, k = 1,2,3, 
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hence = 0 if dimspan{A^^}|^^ = 2, i.e. in the non-degenerate case which is assumed to hold. This 
further yields fjP constant, and there remains the Kirchhoff law 

k=l 

where kP = —divclc = Vct denotes the curvature vector of the contact line C. 

So, Eq = 0 on implies the following: 

(1) The densities are constant, and all velocities vanish. 

(2) The curvatures of the hypersurfaces are constant. 

(3) Ylk=i1^~ , where the coefficients 7 -^ are positive constants. 

But this implies that the classical solution coincides with an equilibrium of the system at any t G 
hence remains at a fixed equilibrium for all t > ti. Consequently, it holds that for any classical 
solution, Eq is strictly decreasing outside of equilibria, i.e. Eq is actually a strict Lyapunov function. 

□ 


To identify all possible equilibrium configurations is a purely geometrical problem. It appears to be a 
challenging problem and will not be analyzed any further, here. 

Final Remarks. 1. The proof that E^ is non-increasing along classical solutions requires all interfacial 
and triple line conditions, in particular the condition ( 66 ). This confirms that the original interface 


formation model of Shikhmurzaev misses one contact line condition. The origin of this transmission 
condition is the fact that transfer of mass, here from one interface across the contact line to another 
interface, is a dissipative process which requires a closure relation. This is similar to the case of mass 
transfer across a fluid interface: even without interfacial mass, a fluid interface carries interfacial energy 
and, in general, entropy can be produced at the interface. In order to avoid entropy production, the 
condition which guarantees zero interfacial entropy production has to be added, leading in the simplest 
case to continuity of the chemical bulk potentials. The triple line analog is equation ( 66 ) above. In 


the more general case of non-trivial entropy production, thermodynamically consistent closure leads 


to a condition like (55). 


2. The molecular kinetic theory of dynamic contact lines supports a friction-like dissipation term at 
the contact line, modeled as being proportional to the square of the contact line speed. If, instead of 


the non-linear closure (55), a linear relation is imposed, the rate of entropy production due to transfer 
of interfacial mass across the contact line becomes proportional to the contact line speed squared. In 
an isothermal setting, this entropy production is proportional to the dissipation of available energy. 
Hence, the so-called contact line friction can be identified with the interfacial mass transfer dissipation 
mechanism. 
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